Pions in the quark matter phase diagram 
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Abstract. The relationship between mesonic correlations and quantum condensates in the 
quark matter phase diagram is explored within a quantum field theoretical approach of the 
Nambu and Jona-Lasinio (NJL) type. Mean-field values in the scalar meson and diquark 
channels are order parameters signalling the occurrence of quark condensates, entailing 
chiral symmetry breaking (xSB) and color superconductivity (2SC) in quark matter. We 
investigate the spectral properties of scalar and pseudoscalar meson excitations in the phase 
diagram in Gaussian approximation and show that outside the x^B region where the pion 
is a zero-width bound state, there are two regions where it can be considered as a quasi- 
bound state with a lifetime exceeding that of a typical heavy-ion collision fireball: (A) 
the high-temperature xSB crossover region at low densities and (B) the high-density color 
superconducting phase at temperatures below 100 MeV. 
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1. INTRODUCTION 

The study of the QCD phase diagram is a key issue in modern theoretical and 
experimental physics of dense matter. Recent heavy-ion collision experiments at 
RHIC Brookhaven [l[ have led to the insight that the quark-gluon plasma (QGP) 
at high temperatures behaves as a perfect fluid with a low viscosity to entropy ratio 
j|/sw0.1-0.2 UMM 

which is very close to the KSS bound [6( for this number, 
l/(47r). This strong deviation from the behavior of a gas of weakly interacting 
quarks and gluons is attributed to the occurence of mesonic bound states @, 0, |jj 
or resonances 0, @, @] in the strongly coupled QGP (sQGP). It has been pointed 
out [To(] that this situation in hot and dense QCD matter bears similarities with 
strongly coupled plasmas in other systems where bound state dissociation or 
Mott-Anderson derealization [ll[ occurs since the effective coupling strength is 
modified by electronic screening and/or Pauli blocking effects. It is thus a very 
general effect expected to occur in a wide variety of dense Fermi systems with 



attractive interactions [TJ, UJ, UJ, [15|, [16|, [17|, 113, 113, 120, [21|. When this Mott 



transition from truly bound to resonantly paired states occurs under conditions 
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of Bose condensation one speaks of a BEC-BCS crossover 22|, |23|, |24(. Recently, 
this transition became accessible to laboratory experiments with ultracold gases 
of fermionic atoms coupled via Feshbach resonances with a strength tunable by 



applying external magnetic fields [2J, [25 



23, 281. The BEC-BCS crossover 



transition in quark matter is of particular theoretical interest due to the additional 
relativistic regime it offers 

Theoretical concepts explaining the appearance of nonperturbative phenomena 
like quantum condensates and bound states in dense Fermi systems with their 
observable consequences shall apply here but must be formulated within a quantum 
field theoretic approach. A systematic treatment of these effects is possible within 
the path integral formulation for finite-temperature quantum field theories. This 
approach is especially suited to take into account the effects of spontaneous 
symmetry breaking. Here we will apply this approach on the example of a model 
field theory of the NJL type to quark matter as a relativistic strongly interacting 
Fermi system. It is our aim to delineate relationships between the regions of the 
NJL model phase diagram where x^B and color superconductivity occur and the 
possibility to observe quasi-bound pionic states in HIC experiments where hot, 
dense QCD matter produced in the form of rather short-lived fireballs. 



2. SCALAR-PSEUDOSCALAR MESONS IN A 
SUPERCONDUCTING TWO-FLAVOR NJL MODEL 

As a generic model system for the description of hot, dense Fermi-systems with 
strong, short-range interactions we consider quark matter described by a model 
Lagrangian with four-fermion coupling. The key quantity is the partition function 
Z from which all thermodynamic quantities can be derived. In the imaginary time 



formalism (t = —it) it can be expressed as a path integral [31 



Z = J V(iq ] )V(q) el" **V-rt<0 , (1) 

where the chemical potential fi is introduced as a Lagrange multiplier for assuring 
conservation of baryon number. The notation d 4 x is shorthand for Jq dr J d 3 x 
where (3 is the inverse temperature. The quark matter is described by a Dirac 
Lagrangian with internal degrees of freedom (Nf = 2 flavors , N c = 3 colors), with 
a current-current-type four-fermion interaction inspired by one-gluon exchange. 
After Fierz transformation of the interaction, we select the scalar diquark channel 
and the scalar-pseudoscalar meson channels so that our model Lagrangian assumes 
the form 

C = q{i$ -m )q + G D {qi^CT 2 \ 2 q T ) {q T iC^T 2 \ 2 q) + G s [ {qqf + {qi^rq) 2 ] ■ (2) 

Here 7^ are the Dirac matrices, A2 is a color SU(3) Gell-Mann matrix, Tj are SU(2) 
flavor matrices and C = ry 2 7° is the charge conjugation matrix. Gs and Go are 
coupling strengths corresponding to the different channels, see for a recent 



review. For the numerical analysis we adopt parameters from Ref. 33] and consider 
Go as a free parameter. 



After introduction of the Hubbard-Stratonovich 34| auxiliary fields A(r, x), 
A*(r,x), tt(t,x), a(r,x) and the Nambu-Gorkov spinors ^ = -^(qq c ) T ', ^ = ^(qq c ) 
with q c (x) = Cq T (x), the partition function becomes a Gaussian path integral in 
the bispinor fields which can be evaluated and yields the fermion determinant 

Z = JvA*VAVaV7re'^ d4xE ^ + ^ -DetlS- 1 ] , (3) 

where the inverse bispinor propagator is a matrix in Nambu-Gorkov-, Dirac-, color- 
and flavor space. After Fourier transformation it reads 

y iA*^r 2 X 2 (iuj n - /j,)>y -m-i~f-p + i>y 5 T t -7T J 

with m = mo + o". So far we could derive with ([3]) a very compact, bosonized form 
of the quark matter partition function ([1]) which is an exact transformation of (pQ), 
now formulated in terms of collective, bosonic fields. As we will demonstrate in 
the following, this form is suitable since it allows to obtain nonperturbative results 
already in the lowest orders with respect to an expansion around the stationary 
values of these fields. In performing this expansion, we may factorize the partition 
function into mean field (MF), Gaussian fluctuation (Gauss) and residual (res) 
contributions 

Z(fi,T) = e-^V = Z MF (fi,T)Z Gauss (v,T)Z rcs (^T) . 

In the following we will discuss the physical content of these approximations. 
In thermodynamical equilibrium, the mean field values satisfy the stationarity 
condition that the thermodynamical potential Qmf = — ^ln^MF be minimal, 

i.e. = = |^ = . This is equivalent to the fulfillment of the gap 

0CTMF OTTMF OAMJf ^ 1 

equations omf = —AGsTt(Smf) = m — m , t^mf — ~ AiGgTr (jstSmf) = and 
Amf = 4Gi)Tr (j5T 2 X 2 Smf) = A, together with the stability criterion that the 
determinant of the curvature matrix formed by the second derivatives is positive. 
After the evaluation of the traces and the sum over the Matsubara frequencies one 
gets 



1 1 -7 (m-m ) 2 |A| 2 1 / . 



pV AGs 4Gd (3V 

(m — m ) 2 |A| 2 f d 3 p 



+ iGb- 4 /(^[^ + + ^ + ^ + 2r " l(1+e " 



AGs AGr> J (2tt) 
+ 2Tln(l + e- /3£ p)+Tln(l + e- /3f p)+Tln(l + e _ ^)l , (5) 



where we have defined the particle dispersion relation E^ = J ' + A 2 with 



£p = E p ± /i, E p = ^m 2 + p 2 . The A^O dispersion law E~ is associated to the 
red and green quarks (antiquarks), whereas the ungapped blue quarks (antiquarks) 
have the dispersion £~ (£+ ). With the aformentioned stationarity conditions applied 




FIGURE 1. Left panel: order parameters for xSB (full lines) and color superconductivity 
(dashed lines) at T = for different values of rjjj. First order phase transitions turn to second 
order or even crossover when rjo is increased. For details, see text. Right panel: phase diagram 
of two-flavor quark matter with critical lines for chiral symmetry breaking (dashed) and color 
superconductivity (solid) for three values of the diquark coupling strength: r/o = 0.75, f .0 and 
1.3. The dots indicate the critical cndpoint for first order phase transitions. 



to Eq. (jSJ) we obtain the gap equations for the order parameters m and A, which 
have to be solved self-consistently, 



m — m 



8G s m 



d 3 p 1 



A 



8G 



D 



+ 

d 3 p 

W) 



l-2n F (E p 



1p_ 

E- 



l-2n F (E-) + l-2n F (E+) 



E: 



E: 



(6) 
(7) 



The Fermi distribution is n F (E) = (1 + e l3E )~ 1 . Solutions of the gap equations for 
the dynamically generated quark mass m and for the diquark pairing gap A at 
T = as a function of fi are shown in the left panel of Fig. [lj From the solutions 
for the order parameters in dependence of the thermodynamical variables T and 
\i we have constructed the phase diagram of the present quark matter model in 
the T — [x plane, see right panel of Fig. [U The two order parameters, that are 
indicators of phase transitions, allow to distinguish 4 phases: A = 0, m ~ mo: 
normal phase (NQM); A^0,m~ mo: color superconductor (2SC); A = 0, m ^> mo: 
chiral symmetry broken phase (%SB); A ^ 0, m > mo: coexistence of %SB and 2SC 
phases. 

Increasing the diquark coupling r/ D = G^/Gs leads to an increase of the diquark 
gap and therefore a rise in the critical temperature for the second order transition 
to a normal quark matter phase. It shifts also the border between color supercon- 
ductivity (2SC) and chiral symmetry broken phase (xSB) to lower values of the 
chemical potential. For very strong coupling t]d ~ 1, a coexistence region developes 



where both order parameters are simultaneously nonvanishing. Since the phase 
border is not of first order, no critical endpoint can be identified in this case. In 
the xSB phase pion and diquark exist as zero width bound states. At the chiral 
symmetry restoration transition, they merge the continuum of unbound states and 
turn into (resonant) scattering states. When this Mott transition occurs within the 
2SC phase we speak of a BEC-BCS crossover. 

Now we discuss the interesting question of the quasiparticle excitations in these 
phases. To this end, we will expand the action functional in the partition function 
up to quadratic order in the mesonic fields and arrive at a tractable approximation 
for the bosonized quark matter model. 

Let us expand now the mesonic fields around their mean field values. In these 
lectures, we will focus on fluctuations in the mesonic channels, where the pion 
and the sigma meson will emerge as quasiparticle degrees of freedom. On the 
example of the pion we will explain the physics of the Mott transition. The detailed 
investigation of the quantized diquark fluctuations, which are also a prerequisite 



of the formation of baryons, will be given elsewhere [3Cj, |35|, |36 . 

As mentioned above, the pion mean field vanishes, but the sigma field can be 
separated into mean field and fluctuations as o — ► a MF + a. Hence it is possible to 
decompose the inverse propagator S~ x into a mean field part and a fluctuation part 
S^ 1 = S MF + T,, so that in the Gaussian approximation the fermion determinant 
becomes 

1 f d 4 q 



Det [S- 1 ] 


Gauss 


Det 


■<y-i " 

°MF 





expj --J j^Tr[SMF(p)Z{q)S M F(p + q)Z{q)] \ , (8) 



(2tt 

where the matrix £ and the propagator Smf are defined as 
_ _ / -a- 275T • 7T 



—O — 275T • 7T 




Smf = I -171+ n- ■ ( 9 ) 



The matrix elements of the Nambu-Gorkov propagator are 

s P tp Po-tp-^P 

where s p ,t p = ±1, (A + ,A~) = (A*, A). For the subsequent evaluation of traces 
in quark-loop diagrams, it is convenient to use this notation with projectors 
in color space, V vg = diag(l,l,0), = diag(0,0,l) and in Dirac space, Ap = 

\ 1 ±70 ^ TP+ m ^ _ xhe summation over Matsubara frequencies po = iuo n is most 
systematic using the above decomposition into simple poles in the p plane. 
The poles of the normal propagators G ± are given by the gapped dispersion 
relations for the paired red-green quarks (antiquarks), = y / (^ p t ) 2 + A 2 , and the 
ungapped dispersions £p = E p ±fj, for the blue quarks (antiquarks). The anomalous 



propagators Fp are only nonvanishing in the 2SC phase where the pair amplitude 
is nonvanishing. Let us notice explicitly that this procedure has yielded an effective 
action that includes the fluctuation terms responsible for the excitation of scalar 
and pseudoscalar mesonic modes. The evaluation of the traces (jHJ) can be performed 
with the result 



-Tr (Smf^Smf^) 



n 



7T,<7 



7V 7T 










(12) 



F~F + + G~G~ +F + F~ I 

x p p+q ~ ^p ^p+q ' p L p+ql ' 



(13) 



where we have introduced the polarization functions 

n CT(T (go,q) = Tr[G+G+ 
n_(g ,q) = -Tr[G p + ( 75 T)G p + +(/ (75T) + F p (75T')F+ +9 (75T) 

+F+( 75 T)F; +g ( 7 5T') + G p (75T*)G; +(? (75r*)] , (14) 



as the key quantities for the investigation of mesonic bound and scattering states 
in quark matter. Indeed by using Bethe Salpeter equation and by evaluating the 
spectral functions one can obtain important information on the mesons properties. 
In the following we perform the further evaluation and discussion for the pionic 
modes, the a mode is treated in an analogous way. We start with the evaluation 
of traces and Matsubara summation and obtain 



n^(g ,q) = 2 



d3 p (i i p (p+q 

2^ \l + SpS k - 



m 



(2tt 



Sp,S k \ 
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Qo + Sfc£p+ q — s p ^p p tp j k E p p Ep +cl q — t k Ep +(l + t p Ep p 
x (tphE^E^ + SpS^^- \Af 
For a pionic mode at rest in the medium (q = 0) this reduces to 



(15) 



lWfc,0) = 8/^{j\T(£,£ 
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(16) 



where we have introduced the phase space occupation factors N(x,y) = l — n F (x) — 
n F (y) (Pauli blocking) and M(x,y) = n F (x) — n F (y). 

We make use of the Dirac identity lim^o^ - = V\ —mb{x) in order to de- 
compose the polarization function into real and imaginary parts after analytical 
continuation to the complex plane. The imaginary part is straightforwardly inte- 
grated after transformation from momentum to energy uj. At the pole, the integra- 
tion variables transform as p u = J < ^~^^^0^ ~ m2 > an d the integration borders 




FIGURE 2. Left panel: Mass spectrum of mesons (tt, a) as a function of the temperature 
for vanishing chemical potential fis — and strong diquark coupling r\u — 1.0. The threshold 
E t h — 2m q for Mott dissociation of pions and occurrence of a nonvanishing decay width T n = 
Im U„/m„ is reached at T Mott = 212.7 MeV (see inset). Right panel: Mass spectrum of mesons 
(tt, a) as a function of the temperature for finite chemical potential fis — 320 MeV and strong 
diquark coupling r/jj — 1.0 in the 2SC phase. Below the threshold E t h = 2 m q for the onset of the 
decay width r 2m there is another process due to the lower threshold E + — E_ switching on (see 
inset). 



shift p G (0,oo) — > uj G (X±,oo), where the thresholds are given by 2m for A = 
and X± = \J(m + fi) 2 + A 2 ± yjm — /i) 2 + A 2 otherwise. In this way one can decom- 
pose the pion polarization function in the 2SC phase into a real and an imaginary 
part. From the real part we can calculate the pion mass by solving the Bethe- 
Salpeter equation 1 — 2G , .Reil 7r7r (m^, 0) = 0, while a nonvanishing imaginary part 
corresponds to a finite pion width for decay into quark-antiquark pairs. 

In Fig. [21 we show the temperature dependence of the masses and widths of 
pions and sigma-mesons for strong diquark coupling t]d = 1.0 at vanishing chemical 
potential (left panel) and at \i B = 320 MeV (right panel). At \i B = 0, the only 
threshold for the imaginary parts of meson decays is 2m, since A = 0. The a mass 
is always above the threshold and therefore this state is unstable in the present 
model. The pion, however, is a bound state until the critical temperature for the 
Mott transition TMott = 212.7 MeV is reached. As can be seen from the slow rise of 
the decay width T n , the pion is still a well-identifyable, long-lived resonance. The 
detailed analytic behavior of the pion at the Mott transition has been discussed in 
the context of the N JL model by Hiifner et al. (3?J , see also the inset of the left panel 
of Fig. [3 It shows strong similarities with the behavior of bound states of fermionic 
atoms in traps when their coupling is tuned by exploiting Feshbach resonances in 
an external magnetic field, see [38( . In the context of RHIC experiments, one has 
discussed such quasi-bound states as an explanation for the perfect liquid behavior 
ofthesQGPfl. 

Next we want to discuss the pionic excitations in the presence of a diquark 
condensate in the 2SC phase, see the right panel of Fig. [21 We observe the 
remarkable fact that the 2SC condensate stabilizes the pion at T = as a true 



H [MeV] 



FIGURE 3. Phase diagram of two- flavor quark matter for the diquark coupling strength 
7/£> = 1.0. The pion is a stable bound state in the horizontally hatched region below the dash- 
dotted line which is correlated with critical line for chiral symmetry breaking (dashed). In the 
vertically hatched region above the dotted line the pion is a short-lived resonance with a width 
rv > 50 MeV < 1 fm/c). There are two regions where the pion is a quasi-bound state in the 
quark plasma: (A) in the crossover region above the chiral phase transition and (B) in the color 
superconducting phase delimited by the solid line. 

bound state, although the pion mass exceeds by far the threshold 2m. This effect 
is due to a compensation of gapped and ungapped quark modes and has been 
discussed before by Ebert et al. [35[ for T = only. Here we extend this study to 
the finite temperature case, where the pion obtains a finite width but is still a very 
good resonance. 

Finally, an interesting insight can come looking at Fig. [3j In this figure we have 
reported the essential phase diagram of the two-flavor quark matter in the case of 
coupling i]£> = 1.00. The dashed line represents the chiral crossover line and the 
solid line refers to the phase transition towards 2SC phase. The dot-dashed line 
shows the border above which the pion turns out to be a quasi bound state in the 
quark plasma. The dotted line indicates the region border (A) below which the 
quasi-stable pionic states have a lifetime greater than 1 fm/c, which is a typical 
lifetime of fireballs @. Thus they would be measured as bound states and could be 
significant in the framework of HIC experiments. Region (B) confirms our claim 
about the absence of stable pions in a two-flavor superconductor. However this 
claim needs a comment. First of all, it can be easily understood how the presence 
of a finite diquark gap (2SC phase) at T=0 stabilizes the pion. Indeed as soon as 
the quark mass drops due to chiral symmetry restoration, the diquark gap tends to 
be finite and thus takes over the role of the quark mass in the dispersion relations. 
Roughly speaking the quark start to be dressed by his interactions. Thus the pion 



2 Note that the physical width T 
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dRcU 
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imn„ 



does not "feel" the drop of the quark mass. But as A melts with increasing T, the 
pion width raises, first slowly and then rapidly, reflecting the behaviour of A as a 
function of T. This leads inevitably to a destabilization of the the pion states in the 
vicinity of the 2SC phase border. The discussion of the mesonic modes in the 2SC 
phase points to a very rich spectrum of excitations which eventually leads to specific 
new observable signals of this hypothetical phase. The CBM experiment planned 
at FAIR Darmstadt and the NICA project at JINR Dubna could be capable of 
creating thermodynamical conditions for the observation of these excitations in 
the experiment. One promising signal could be the scalar resonance in the pion- 
pion scattering at the two-pion threshold which is in principle observable, e.g., in 



the two-photon decay channel |39l . |40| . However, the description of this state goes 



beyond the Gaussian approximation to which we restrict ourselves in this work. 



3. CONCLUSIONS 

In this work we have derived and evaluated the gap equations and the scalar- 
pseudoscalar meson spectra within a path integral approach to the two-flavor NJL 
type model of superconducting quark matter. 

After fixing the parameters of the model to the light meson spectrum in the 
vacuum, the diquark coupling remains as a free parameter which has been used to 
extend the model beyond the traditional range of applications into the region of 
BEC-BCS crossover. 

We have presented the phase diagram of quark matter at strong and very strong 
coupling. The origin of the BEC-BCS crossover in superconducting quark matter 
is the Mott transition for diquark bound states. We explain the physics of the Mott 
transition on the example of mesonic correlations. We have investigated the meson 
spectra (bound and scattering states) outside (T > T c ) and for the first time also 
inside (T < T c ) the color superconductivity region. We find the thresholds for the 
dissociation of pionic bound states into unbound, but resonant scattering states in 
the quark-antiquark continuum and we have shown that outside the %SB region 
where the pion is a zero-width bound state, there are two regions where it can 
be considered as a quasi-bound state with a lifetime exceeding that of a typical 
heavy-ion collision fireball: (A) the high-temperature %SB crossover region at low 
densities and (B) the high-density color superconducting phase at temperatures 
below 100 MeV. 
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